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We study the localization transitions which arise in both one and two dimensions when quan- 
tum mechanical particles described by a random Schrodinger equation are subjected to a constant 
imaginary vector potential. A path-integral formulation relates the transition to flux lines depinned 
from columnar defects by a transverse magnetic field in superconductors. The theory predicts that 
the transverse Meissner effect is accompanied by stretched exponential relaxation of the field into 
the bulk and a diverging penetration depth at the transition. 



PACS numbers: 05.30.Jp, 72.f5.Rn, 74.60. Ge 

Although forbidden in conventional quantum mechan- 
ics, exponentiated non-Hermitian quantum Hamiltonians 
do appear in the transfer matrices of classical statistical 
mechanics problems. A nonequilibrium process can be 
described as the time-evolution of a non-Hermitian sys- 
tem [[[] . Another example is the XX Z spin chain mapped 
onto the asymmetric six- vertex model . 

In this Letter, we investigate a non-Hermitian quan- 
tum Hamiltonian with randomness. The study is moti- 
vated by a mapping of flux lines in a d+ 1-dimcnsional su- 
perconductor to the world lines of ci-dimensional bosons. 
Columnar defects in the superconductor, which were in- 
troduced experimentally in order to pin the flux lines || , 
give rise to random potential in the boson system Q . Al- 
though the field component H z parallel to the columns 
acts like a chemical potential for the bosons, the compo- 
nent perpendicular to the columns results in a constant 
imaginary vector potential [^) . 

We study localization in this simple example of non- 
Hermitian quantum mechanics, and thereby show how 
a flux line is depinned from columnar defects by an in- 
creasing perpendicular magnetic field H±. It is gener- 
ally believed that all eigenstates are localized in conven- 
tional one- and two-dimensional non-interacting quan- 
tum systems with randomness. On the other hand, it 
is almost obvious that a flux line is depinned from de- 
fects by a strong perpendicular field component. This 
indicates that the present non-Hermitian system has ex- 
tended states in a large H ± region, and that there must 
be a derealization transition at a certain strength of H j_ . 

The non-Hermitian Hamiltonian treated hereafter has 
the form H = (p + ih) 2 /(2m) + V(x), where p = 
(h/i)W, and V(x) is a random potential. The non- 
Hermitian field h originates in the transverse magnetic 
field as h = 4>qH j_/(47t), where 0o is the flux quantum 
HI). Figure 1 shows a vortex whose "world line" is de- 
scribed by this Hamiltonian with periodic boundary con- 
ditions in one dimension. The mass m is equivalent to 
the tilt modulus of the flux line. The Planck parameter 
h corresponds to the temperature of the superconduc- 
tor, while the inverse temperature of the quantum system 
corresponds to the thickness L T of the superconductor. 
Interactions between many particles (or flux lines) can be 



treated approximately by forbidding multiple occupancy 
of localized state in a tight binding model (see below) [Q . 
Interactions in the delocalized regime will be discussed in 
a future publication Q. 

Since the field h acts like a vector potential, we define 
the current operator as J = —idH/dh = (p + ih)/m. 
The imaginary part of the current describes the tilt slope 
of a flux line. To see this, note first that the posi- 
tion of the flux line x at the distance r from the bot- 
tom surface of the superconductor is given by (x) T = 
Z-^^le-^-^^xe-^^l^), where V* and ^ de- 
scribe boundary conditions at the bottom and top sur- 
faces (r = 0, L T ) of the superconductor, respectively. 
The partition function is Z = (<0 / |e~ LT ' H / fi |^ i ). The 
commutation relation [H, x] ~ ~ihJ leads immediately 
to (d/dr) (x) T = -i{J) T =Im{J) T . 

The total displacement of the flux line between the 
bottom and the top surfaces is given by (x) L — (x) = 

h(d I dh) In Z = Ira fy T (J) t cIt. This quantity is an in- 
dicator of the delocalization transition; the transverse 
displacement of a flux line must be order of the system 
size when it is depinned. 

Let us first consider localized states in a small perpen- 
dicular field. Assume the eigenfunctions ipn( x ) and the 
eigenvalues e n are known for h = 0. For small h the 
right- and left-eigenvectors of Ti. are given by 4>ni x ', h) = 
e h - x / h 4> n {x-h = 0), and ip%(x;h) = e^^ip^x; h = 
0). The energy eigenvalue e„ is unchanged under this 
"imaginary" gauge transformation Jt|. The imaginary 
gauge transformation applies even to many-body eigen- 
vectors of interacting systems (6). However, the above 
wave functions ipn and ma y diverge in \x\ — > oo, and 
hence may not be normalizable. The normalizablity con- 
dition is | ft. | < TiK n , where n n is the inverse localization 
length of the state tp n (x;h = 0). Then the wave func- 
tion (with the normalization J d d xtp^(x)ip^ (x) = 1) is 
approximately 

where x n is the localization center for h — and fid is 
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FIG. 1. One flux line (wavy curve) induced by the field H z 
and interacting with columnar pins in a cylindrical super- 
conducting shell with radial thickness smaller than the pen- 
etration depth of the defect-free material. The field H± is 
generated by the current I threading the ring. 

the total solid angle of the d-dimensional space. 

We naturally regard the point \h\ — hn n as the der- 
ealization point of the state ip^- ^ n the region \h\ > hn n , 
the eigenfunction becomes extended. Hence we need to 
specify the boundary conditions in order to obtain a well- 
defined wave function in the thermodynamic limit. 

Let us consider imposing the periodic boundary con- 
dition %f)^(L x /2,y, . ..) = ip*(-L x /2,y, . . .), with the x 
axis parallel to h. The one-dimensional periodic system 
is realized in the setup shown in Fig. 1. Except for the 
imaginary vector potential, this setup is the flux-line ana- 
logue of a mesoscopic metal ring threaded by a solenoid. 
At the boundaries x — ±L x /2, the wave function of 
the form ([!]) has a mismatch of the factor e ^^~ h l n ) L ^ . 
In the region h < hn n , this mismatch is exponentially 
small, and so is the change necessary to meet the periodic 
boundary condition. In the region h > hn n , however, the 
wave function changes drastically, and a complex eigen- 
value appears. To understand this point, consider the 
case \h\ — > oo, a limit in which the random potential 
V(x) may be neglected [||. The periodic boundary con- 
dition is satisfied only by the extended function e l ^' x 
with k v = 2n v ii / L u , where n v is an integer and L v is the 
system size in the x v direction. (If a right-eigenvector is 
e ik x^ ^he correS p 0nc ii n g left-eigenvector is e~ l ^' x .) The 
eigenvalue is complex: e(k) — (hk + ih) 2 /(2m). 

The imaginary part of an eigenvalue has the follow- 
ing meaning. The time evolution of an extended state 
is described as e ^k-x-re/K _ This must actually be a 
function of k ■ (x — rim (J)), because the correspond- 



ing flux line has the tilt slope Im(J). (The angular 
brackets here denote the expectation value with respect 
to the state with energy e.) We hence expect the relation 
Ime ~ hk ■ Im ( J) = — hk ■ (d/dh)Ree. Thus appear- 
ance of a complex eigenvalue is another indicator of the 
delocalization transition. 

It is instructive to solve the one-dimensional sys- 
tem with a single attractive impurity: TL = — (KV — 
h) 2 /(2m) - V 6(x). We first solve H<f>(x) = e<j)(x) for 
x 7^ to obtain the two degenerate solutions <f>±(x) = 
e ±ikx+hx/h ^ anc j ^en write the general solution as ip = 
A(f> + + B<f>- . We impose the periodic boundary condi- 
tion ip(L x ) — ip(0) and the condition for the delta po- 
tential, tp'(L x ) - f(0) = (2mV /h 2 )xj;(Q). A nontriv- 
ial solution arises whenever k [cosh(L x h/K) — cos(L x k)] + 
(mVo/h 2 ) sin(L x k) = 0. The ground state (defined by 
the lowest real part of the eigenvalue) is localized only 
for h < hK gs — mVo/h. All the other states are ex- 
tended. The localized ground-state energy is e gs = 
-itiVq /(2h 2 ), and the wave function takes the form 
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tion length in the direction of h grows as (hn gs — h) , 
and the state goes through a delocalization transition at 
h = hn gs . The ground state for h > hK gs is extended 
with the energy —h 2 /(2m). The corresponding flux line 
has the tilt slope h/m, which is equal to the value in the 
impurity-free case. The tilt slope of the ground state has 
a jump at the delocalization point. 

Extended states in the thermodynamic limit have 
a form close to the pure case. The leading term is 
" aX /\/Lx with k n — nir/L x , where n is odd for h < 

They also have the re- 
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gs , and is even for h > hn gs . 



flection term e -i- k n x + 2hx [ n the region x < 0. The energy 
eigenvalue takes the same form as the pure system. The 
above solution shows that a depinned flux line almost 
ignores the columnar defect. 

Now we move on to the case of random potential. For 
the purpose of numerical calculations, it is convenient to 
introduce the non-Hcrmitian tight-binding model. The 
second-quantized Hamiltonian is written in the form 



H 



x 



X v=l 

VxbLbx, 



x+e v u x 



e "x^x+e, 



(2) 



where the {b x ,bx} are boson creation and annihi- 
lation operators g, and the {e„} are the unit lat- 
tice vectors. The hopping element is approximately 
given by g] t ~ V hind ex.p(-^/2mV hind a/h), where 
Vbind is the binding energy in the tight-binding ap- 
proximation, and a is the lattice spacing. We apply 
the periodic boundary conditions bx+N v e v = bx for 
v = 1,2, ... ,d, where N v = L v /a. The complex eigen- 
values in this non-Hermitian system always appear in 
conjugate pairs; if there is a complex eigenvalue e with 
a right-eigenvector tp R , there is also the eigenvalue e* 
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FIG. 2. (a) Ime vs. Ree for a particular sample of the d = 1 
tight-binding model with A/t = 1 and = 1000a. Each 
eigenstate is marked by a cross +. Plots for different h are 
offset for clarity, (b) Im J (fluxon tilt) vs. Re e for the same 
sample as in (a). The inset shows na vs. Ree for A/t = 1, 
L x = 500a and h — 0, averaged over one hundred samples 
binned with energy window 0.04i. 

with the right-eigenvector (if; R )* . This symmetry ensures 
that the partition function 2 is real. Another symmetry 
is Ti(h) T = H(—h). Because of this symmetry, a right- 
eigenfunction of TL(h) equals to the left-eigenfunction 
of H(—h) with the same eigenvalue. If Vx = 0, the 

eigenvalues are e — — * X^=i cos [(^ + ih u /K)a] with 
k v = 2n„7r/L„, where n v is an integer. The eigenfunc- 
tion takes the form ip R (x) oc exp(ifc ■ x). Numerical 
calculations for Vx ^ were carried out assuming that 
the random potential Vx is uncorrelated in space, and 
uniformly distributed in the range [—A, A]. Results for 



random hopping models, with no site randomness, are 
qualitatively similar || . 

Figure 2 (a) shows the eigenvalue spectrum for d = 1. 
The real eigenvalues indicate localized states, while the 
complex eigenvalues indicate delocalized ones. As dis- 
cussed above, the eigenvalues for the localized states are 
indeed independent of h. The behavior of the delocalized 
states, on the other hand, is similar to impurity-free sys- 
tems except near the mobility edges. Figure 2 (b) shows 
the imaginary part of the current, or the tilt slope of the 
corresponding flux line plotted against the real part of 
the energy. The tilt slope has a jump at the mobility 
edges as in the one-impurity case. The inset to Fig. 2 
(b) shows the inverse localization length for this problem 
when h — 0. According to the derealization criterion 
h c = Kk, a region of delocalized states should appear 
as h increases first at the band center and then move 
outwards, as we observe. 

In two dimensions we again find bands of localized en- 
ergies bounded by a mobility edge, although extended 
and localized states are apparently mixed in a compli- 
cated way near the band center. Extended states again 
almost ignore impurities. The case of an attractive im- 
purity shows that the localized ground state exists in a 
certain region of small h, and that the energy spectrum 
of the extended states is close to the impurity-free case. 

An interesting question in the two-dimensional random 
case is whether a delocalized state is extended in both the 
x and y directions, or extended only in the direction of 
h. We conclude on the basis of the following argument 
that the former is the case. If the field h is parallel to 
the x axis and a state is extended in this direction, we 
may assume that the wave function has the approximate 
form L x 1 ^ 2 e lk * x (f>(y) in order to accommodate periodic 
boundary conditions. This is indeed an eigenfunction of 
the partial Hamiltonian which has only the kinetic term 
in the x direction. We then calculate the effective Hamil- 
tonian for <fi(y) in degenerate perturbation theory, taking 
the kinetic term in the y direction and the random poten- 
tial term as the perturbation. We find that the effective 
Hamiltonian for <fr(y) has a random potential term of the 
form L~ x J V(x,y)dx. The width of the probability dis- 
tribution of this effective random potential vanishes as 

— 1/2 

L x in the limit L x — > oo. The flux line is a random 
walker in the y direction, and hence <fi(y) is an extended 
state. Derealization in directions both parallel and per- 
pendicular to h is more readily observed numerically for 
h along the diagonal of a square lattice than for h parallel 
to e x or e y [§ . 

Spectra such as those in Fig. 2 (a) at intermediate val- 
ues of h are relevant at low but finite concentrations of 
interacting flux lines, provided that we fill up the local- 
ized states in order of increasing energy up to energy 
e = (i, thus forbidding multiple occupancy of a single lo- 
calized state [[| . The system is then characterized by an 
average chemical potential /i = ^(H z ) which controls the 
flux-line density and separates occupied from unoccupied 
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levels. Once this chemical potential exceeds the mobility 
edge, however, we expect spiral trajectories and "bose 
condensation" of delocalized flux lines into the lowest- 
energy extended state; more sophisticated methods of 
dealing with interactions are then required |^j^J . 

Consider the probability distribution of a flux line near 
a free surface. At the distance r from the bottom surface 
of the superconductor this probability may be written 
P(x;t) = Z~ 1 {^\e-( L -~ T ') H l h \x){x\e- Tn l h \ij) i ), where 
we assume free boundary conditions at the bottom and 
the top surfaces: = \ip*) — J d d x'\x'). In the 

limit L T — > oo, the probability distribution of the most 
weakly bound flux line at the top surface is propor- 
tional to the right-eigenvector with the eigenvalue e = (i, 
P(x;L T ) oc {x\ip^) — tp^(x), while the distribution at 
the bottom surface is proportional to the left-eigenvector 
of the same state, P(x; 0) oc (ipn\x) = tpj^(x). The distri- 
bution far from these surfaces is given by P{x; L T /2) — 
ip^(x)ipL(x). In the localized regime, the imaginary 
gauge transformation for and tpj^ describes the dis- 
placement of the flux line at the surfaces, while the bulk 
distribution P(x; L T /2) is independent of h, consistent 
with the transverse Meissner effect 0. These state- 
ments are approximate because we have neglected exci- 
tations of occupied states. They are exact for the ground 
state. The displacement of the flux lines near the top and 
bottom surfaces allows partial penetration of the trans- 
verse component of the magnetic field, characterized by 
a transverse London penetration depth t* . 

In the case of one columnar defect (or a single impurity 
in a tight binding model) , the displacement of a flux line 
takes the form ({x) T — (x)^) ■ h ~ e -rAei/n^ w ] lere 
is the energy gap between the localized ground state and 
the extended first excited state. In the random case with 
finite concentration of flux lines, we consider the relax- 
ation near a free surface r = of the most weakly bound 
line to the highest-energy occupied state. We show below 
that the displacement has the stretched exponential form 

«*) T - (x)^) • h ~ exph a (r/T*) 1/(d+1) ], (3) 

r — >oo 

where (x)oo is the center of the localized state in 
the bulk, the penetration depth is given by t* = 
g(^)h d+1 1 '(Uk(ix) — h) d , and a is a constant of order of 
unity. Here g(/J.) and are respectively the density of 
states and the inverse localization length at the chemical 
potential. Note that r* ~ where £j_ ~ l/(hn — h) is 
the diverging surface localization length near the transi- 
tion. 

For the derivation of the above, we expand P(x;t) 
with respect to the eigenfunctions. In the limit L T — > oo, 
we have P(x;r) ~ (tp fl \x)J2' n c n (x\ipn)e~ TA£ " /h , where 
Ae„ = e„ — //, (ipn\x) = ip^{x) is the left-eigenvector 
describing the shifted probability distribution of the sur- 
face, and the {c„} are certain constants coming from the 
normalization factor. The summation Yl'n 1S restricted 



to states with energies e > /i. The asymptotic form 
(|l|) for localized states now leads to \(x) T — (x)^] ~ 
J2' n c n e~ fn/h , where /„(r) = (hn n - hcos6 n )r n + rAe„, 
r n = \x n - a; M |, cos6»„ = h ■ (x^ - x n )/\h\r n , and 
the {c„} are additional constant factors. Upon setting 
g(n)r d A.£ n ~ 1 |Q, we find that the two terms in /„ 
compete with each other; in other words, the further the 
flux line hops, the more the kink energy costs, but the 
lower the binding energy of the state n is. Minimizing 
/„ with respect to r„ and 9 n , we arrive at Eq. (||). We 
expect that a similar relaxation process governs lines re- 
laxing from the surface to occupied states with e < \l as 
t — ► oo . For d = 1 and randomly distributed column po- 
sitions, rare events like exceptionally small hopping ma- 
trix elements will alter Eq. (|j|) for very large r, similar 
to flux lines pinned by parallel, randomly placed grain 
boundaries [ |10| . 

The behavior of the flux-line tilt modulus (inverse bo- 
son superfluid density) and flux-flow resistivity above the 
mobility edge may be obtained || using the methods of 
Ref. §. 
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